clear all
clc
ft = @(x) 1/sqrt(1+x.^3);
I= mymulNewtonCotes(ft,0,1,10,2)
